Computer Aritmetiek 

Pascal Vyncke 


OPGAVE 

(1) Om 7 beduidende cijfers voor in 2 uit te rekenen door gebruik te maken van de 
reeksontwikkeling van ln(1+x) benodigt men een ongelooijk groot aantal termen. 
Hoeveel termen heb je nodig als je In 2 wil programmeren door gebruik te maken 
van de reeksontwikkeling van In _i+x 

inx_(het is voor het project voldoende om alleen 

de truncatiefout in rekening te brengen)? 

(a) Is deze x exact voorstelbaar? 

(b) Programmeer de berekening van In 2 zo stabiel mogelijk in enkele en dubbele 
precisie. 

(c) Heb je in enkele precisie e_ectief 7 beduidende cijfers of wordt de afrondingsfout 
te groot? 

(2) Voor xo > Dl geldt dat de rij fxngn 2 N met 
Xn+1 = 2n+1 _p1 + 2ünXn L 1_ 
convergeert naar ln(1 + xo). 

(a) Programmeer deze recursie en vergelijk met ln(1 + xo). 

(b) Herschrijf de recursie zodat de berekeningen betrouwbaarder worden. 

(3) We willen f(x) = (x □ sin x)=x3 berekenen voor x = 1Onk met k = 0; : : : ; 15. 

(a) Welke waarden bekom je als je f(x) programmeert zoals opgegeven? 

(b) Waaraan is limxiof(x) gelijk? 

(c) Hoe moet je de berekening van f(x) voor x klein organizeren? 


Oefening 1 

1. a) 

Om de uitkomst van ln2 te verkrijgen via de reeksontwikkeling van in [(l+x)/(l-x)j is x niet 
exact voorstelbaar (niet in binaire representatie, niet in decimale representatie). 

1. b) 

De berekening van ln2 via de reeksontwikkeling van ln [(l+x)/(l-x)j is het volgende. 

Volgens de definitie van ln is het volgende waar: 
ln [(l+x)/(l-x)j = ln (1+x) - ln (1-x) 

We kunnen ln(l+x) en ln(l-x) omzetten in elk hun eigen reeksontwikkeling, en vervolgens 
van elkaar aftre kk en en zo een totale reeksontwikkeling verkrijgen die ln [(l+x)/(l-x)j 
voorstelt. 

Hiermee komen we echter niet verder, aangezien we dan twee keer de ln moeten berekenen en 
aftrekken (en aftrekking is “gevaarlijk” voor computers op basis van nauwkeurigheid; om aan 
de opgave te voldoen, stabiel programmeren dienen we dus de aftrekking te vermijden). 
Bovendien wordt in de opgave ook aan ge geven dat de reeksontwikkeling van ln(l+x) al zeer 
groot is en dat we het hiermee niet gaan halen. 


Dus ik zoek het via een andere weg de reeksontwikkeling van ln [(l+x)/(l-x)j. 




De reeksontwikkeling van ln [(l+x)/(l-x)] is: 

We geven er een functie f(x) als naam aan 
f(x) = ln[(l+x)/(l-x)] 

via afleiden van f(x) krijgen we: 
f (x) = 2 / (1-x 2 ) 

de reeksontwikkeling hiervan is: 

2 + 2x 2 + 2x 4 + 2x 6 + 2x 8 + 2x 10 +... 

door terug te integreren term per term van bovenstaande reeksontwikkeling, krijgen we: 
2x+(2x 3 )/3+(2x 5 )/5+(2x 7 )/7+. ...) 

= 2(x+x 3 /3+x 5 /5+x 7 /7+....) en dit is de reeksontwikkeling van ln [(l+x)/(l-x)] 

We dienen dus 2(x+x 3 /3+x 5 /5+x 7 /7+....) te programmeren en berekenen. 

Deze formule is echter geldig voor eender welke ln. De opgave vraagt enkel de ln 2. 

Omdat de formule enkel geldig is voor een willekeurige x tussen -1 en 1, gaan we dit 
omzetten. Het is namelijk mogelijk om elk getal hierin om te zetten. 

Het getal 2 ligt niet tussen de grenzen van -1 en 1. Echter, door dit te inverteren, dus 2" 1 , 
komt dit wél te liggen in -1 en 1, namelijk Vi. 

We kunnen dit gebruiken omwille van volgende regel van logaritmen: 

-ln(x) = ln(l/x) 

Zodus, kunnen we de ln(2) gaan voorstellen, door gebruik te maken van Vi en dan onze 
uitkomst uiteindelijk het teken om te draaien (dus positief wordt negatief, en omgekeerd). 

Voor gebruik te maken van onze formule, bekomen we dus dat x = -1/3 

Dit getal is (zoals in onderdeel a van deze oefening gezegd), niet exact voorstelbaar. 

(Namelijk in decimale getallen is dit -0,3333333333333333... en ook in binaire voorstelling, 
namelijk-.01010101010101010...) 

We kunnen dus het volgende stellen: 

ln 2 = ln [(l+(-1/3))/( 1 -(-1/3))] = 2((-l/3)+ (-l/3) 3 /3+(-l/3) 5 /5+(-l/3) 7 /7+....) 

Een narekening van dit resultaat in een wiskundig programma (Mathlab), geeft aan dat 
2((-l/3)+ (-l/3) A 3/3+(-l/3) A 5/5+(-l/3) A 7/7) ~ -0,693134757 

De vraag wat ln2 nu eigenlijk is, krijgen we als antwoord van het programma 
0,6931471805599450, waarbij de cijfertjes in het vet dezelfde zijn als met ons vorige resultaat 
verkregen, met andere woorden zitten we daarmee al aan 5 juiste cijfers. 


Hiermee is de oplossing gecontroleerd, en kunnen we nu deze oplossing gaan implementeren. 
Ik heb gekozen om dit in C++ te doen. 



Het stabiel programmeren van de ln2 is reeds gedeeltelijk opgelost. Ik gebruik nu al de 
reeksontwikkeling en er worden enkel maar optellingen gebruikt, en geen aftrekkingen (wat 
bij de allereerste oplossing, namelijk ln [(l+x)/(l-x)] = ln (1+x) - ln (1-x) niet zo was). 

Een tweede regel dat ik hanteer is dat we éérst de kleinste getallen moeten optellen, en dan 
pas de grote. Dus de reeksontwikkeling in feite omdraaien (wat met optellen volgens de 
wiskunde perfect mag). Dus: 

ln [(1+(-1/3))/( 1 -(-1/3))] = 2(....+(-1/3) 7 /7+(-1/3) 5 /5+(-1/3) 3 /3+(-1/3)) 

Dit blijft dus hetzelfde als resultaat hebben wiskundig gezien, met de computer uitgerekend 
verhoogt dit echter de nauwkeurigheid (beter gezegd: laat de (kans op) afrondingsfouten 
dalen). 

We willen 7 beduidende cijfers hebben, we nemen daarom 8 termen. 


De C++ code: 

Voor single precision: 

#ifndef INC_IOSTREAM 
#include <iostream> 

#define INC_IOSTREAM 
#endif 

#ifndef INC_STRING 
#define INC_STRING 
#include <string> 

#endif 

#ifndef INC_MATH 
#define INCJYLATH 
#include <math.h> 

#endif 

using namespace std; 

int main () { 

float x=0; 
float oplossing=0; 
float tussenoplossing=0 ; 
int macht; 

macht = 13; // Komt overeen met 6 termen ((6*2) + 1) 
x = 1; 

x = -(x / 3); // x = -0.33333333333 

while (macht >0) { 

// Een tussenoplossing berekenen 


tussenoplossing = (pow(x,macht) / (macht)); 

// Optellen met onze vorige oplossing 
oplossing = oplossing + tussenoplossing; 

macht = macht - 2; 

} 

oplossing = (oplossing * 2); // Oplossing vermenigvuldigen met 2 

oplossing = -oplossing; // Positief wordt negatief, en omgekeerd 

// Nu onze gevonden oplossing outputten 
cout.precision(50); 
cout << oplossing; 

} 

Dit geeft uitgevoerd op een Pentium 4 met het Microsoft Windows platform het volgende als 
resultaat: 0.69314718246459961 
Hierbij zijn er 7 beduidende cijfers: 

0.6931471 8246459961 


Voor doublé precision 

#ifndef INC_IOSTREAM 
#include <iostream> 

#define INC_IOSTREAM 
#endif 

#ifndef INC_STRING 
#define INC_STRING 
#include <string> 

#endif 

#ifndef INCJYLATH 
#define INCJYLATH 
#include <math.h> 

#endif 

using namespace std; 

int main () { 

doublé x=0; 
doublé oplossing=0; 
doublé tussenoplossing=0; 
int macht; 

macht = 15; // Komt overeen met 7 termen ((7*2) + 1) 
x = 1; 

x = -(x / 3); // x = -0.33333333333 

while (macht >0) { 

// Een tussenoplossing berekenen 
tussenoplossing = (pow(x,macht) / (macht)); 


// Optellen met onze vorige oplossing 



oplossing = oplossing + tussenoplossing; 


macht = macht - 2; 

} 

oplossing = (oplossing * 2); // Oplossing vermenigvuldigen met 2 

oplossing = -oplossing; // Positief wordt negatief, en omgekeerd 

// Nu onze gevonden oplossing outputten 
cout.precision(50); 
cout << oplossing; 

} 

Dit geeft uitgevoerd op een Pentium 4 met het Microsoft Windows platform het volgende als 
resultaat: 0.693147179548241230 
Hierbij zijn er 7 beduidende cijfers: 

0.6931471 79548241230 


l.c) 

In single precision hebben we net voldoende precisie om 7 beduidende cijfers te hebben. 
We kunnen dit uitrekenen met de formule. Namelijk: 


7 beduidende cijfers O 


x 2 


Dit uitgerekend brengt dat: 

7 beduidende cijfers O 0,00000000274797869995 < 0,0000005 
En dit klopt want dit is inderdaad kleiner. 


Dit valt ook te vergelijken: 

Single precision: 0.6931471 82464699610 
Doublé precision: 0.6931471 79548241230 

Juist na wat vetjes gedrukt is, zijn dit andere cijfers. De cijfers van doublé precision zijn 
echter wel correct, die van single niet meer. 


Samengevat: het lukt juist om 7 beduidende cijfers te hebben met single precision. 








Oefening 2 


2a) 

De recursie heb ik geprogrammeerd in C++. 

De C++ programmacode: 

#ifndef INC_IOSTREAM 
#include <iostream> 

#define INC_IOSTREAM 
#endif 

#ifndef INC_STRING 
#define INC_STRING 
#include <string> 

#endif 

#ifndef INC_MATH 
#define INCJYLATH 
#include <math.h> 

#endif 

using namespace std; 

doublé Oefening2( int n, doublé xO) { 
if (n > 0) { 

return ( pow(2,n) * (sqrt(l + ((1 / pow(2,n-l)) * 
Oefening2((n-1), xO)))-1)) ; 

} 

else 

{ 

return xO; 


int main () { 

doublé xO; 


// Bereken ln(1+2) 
xO = 2; 

cout.precision(50); 


cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(10,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(15,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(20,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(25,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(30,xO) 

« 

endl 


cout << endl << endl; 

cout << "Vergelijking met log(l+x0)" << log(l+x0) << endl; 


// Bereken ln(l+3) 
xO = 3; 

cout.precision(50); 


cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(10,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(15,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(20,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(25,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2(30,xO) 

« 

endl 


cout << endl << endl; 

cout << "Vergelijking met log(l+x0)" << log(l+x0) << endl; 

} 
















Bovenstaande programma gaat de ln (1+2) en van ln (1+3) uitrekenen. 


Bovenstaande programma geeft op een Pentium 4 met het Microsoft Windows platform het 
volgende als resultaat: 

Uitkomst voor x0=2 :1.0992018300289601000000000000 
Uitkomst voor x0=2: 1.0986307054554345000000000000 
Uitkomst voor x0=2: 1.0986128642689437000000000000 
Uitkomst voor x0=2: 1.0986123085021973000000000000 
Uitkomst voor x0=2: 1.0986120700836182000000000000 
Vergelijking met log(l+x0) = 1.098612288668109800000000000000 

Uitkomst voor x0=3: 1.3872331695188223000000000000 
Uitkomst voor x0=3: 1.3863236860561301000000000000 
Uitkomst voor x0=3: 1.3862952776253223000000000000 
Uitkomst voor x0=3: 1.3862943947315216000000000000 
Uitkomst voor x0=3: 1.3862943649291992000000000000 
Vergelijking met log(l+x0) = 1.386294361119890600000000000000 


Het eerste wat we opmerken is dat naarmate we meer recursies doen, het resultaat 
nauwkeuriger en nauwkeuriger wordt. In mijn programma ga ik tot 30 stappen. Hiermee kan 
ik, in vergelijking met de “echte” uitkomst van het logaritme, tot 8 cijfers nauwkeurig zijn: 

Uitkomst voor x0=3 na 30 recursies: 1.38629436 49291992 
Vergelijking met log(l+x0) = 1.38629436 1 1198906 


2b) 

In oefening a is de berekening met de recursie puur ingegeven zoals in de opgave, zonder 
rekening te houden met bepaalde optimalisaties, regels of om de nauwkeurigheid te verhogen 
(en afrondingsfouten te vermijden). 

Dit wordt nu wel toegepast, om zo een mogelijk beter resultaat te verkrijgen. 

Origineel volgens opgave: 

=2”'(Vl + 2-*„-l) 

Om dit effectief te programmeren, werd dit omgezet naar hetvolgende: 

Of met andere woorden werd n+1 gelijk aan n. We roepen namelijk in ons programma x n op 
en niet x n +i. 



Met de eigenschappen van de wortels, kunnen we dit omvormen: 



En dan gebruiken we een eigenschap van de wortels om dit te herschrijven: 



We gaan nu de 1 - 1 schrappen in de teller: 



Nu kunnen we de macht 2" (n_1) uit onze vergelijking halen: 



Of met andere woorden: 



En dit is gelijk aan: 



Hierdoor hebben we deze grote machten en delingen weggewerkt, waardoor er een hogere 
nauwkeurigheid wordt verkregen. Dit bovenstaande kunnen we nu gaan programmeren. 



In C++ wordt dit: 


#ifndef INC_IOSTREAM 
#include <iostream> 

#define INC_IOSTREAM 
#endif 

#ifndef INC_STRING 
#define INC_STRING 
#include <string> 

#endif 

#ifndef INC_MATH 
#define INC_MATH 
#include <math.h> 

#endif 

using namespace std; 

doublé Oefening2b( int n, doublé xO) { 

doublé tempO, templ, temp2, temp3, temp4, temp5; 
if (n > 0) { 

return ((((Oefening2((n-1), xO)))/(sqrt(1 + ((1 / pow(2,n-l)) 
Oefening2((n-1), xO))) + 1))*2) ; 

} 

else 

{ 

return xO; 

} 


int main () { 

doublé xO; 


// Bereken ln(1+2) 
xO = 2; 

cout.precision(50); 


cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(10,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(15,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(20,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(25,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(30,xO) 

« 

endl ; 


cout << endl << endl; 

cout << "Vergelijking met log(l+x0)" << log(l+x0) << endl; 


// Bereken ln(l+3) 
xO = 3; 

cout.precision(50); 


cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(10,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(15,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(20,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(25,xO) 

« 

endl 

cout 

« 

"Uitkomst 

voor 

x0 = " 

« 

xO 

« " 

" « 

Oefening2b(30,xO) 

« 

endl ; 


cout << endl << endl; 

cout << "Vergelijking met log(l+x0)" << log(l+x0) << endl; 
















De bekomen recursie, dus: 


x„ = 2. 


V^-'-'v.+i 


Kan terug herschreven worden in de vorm zoals in de opgave: 


X n+1 = 2 -l 


Jl + 2 n x n +1 


Deze geeft betere waarden terug. 



Oefening 3 

3a) 

De te programmeren functie is: 


ƒ(*) = 


x - sin x 


Ook deze oefening maakte ik in C++. De programmacode is als volgt: 

#ifndef INC_IOSTREAM 
#include <iostream> 

#define INC_IOSTREAM 
#endif 

#ifndef INC_STRING 
#define INC_STRING 
#include <string> 

#endif 

#ifndef INC_MATH 
#define INC_MATH 
#include <math.h> 

#endif 

using namespace std; 

doublé 0efening3a (doublé x) { 

return ((x - sin(x))/pow(x, 3)); 

} 

int main () { 

int k = 0; 
doublé x; 

cout.precision (50) ; 

while (k <= 15) { 

x = (1 / pow(10, k)) ; 

cout << "f(x) met x=10 A (-" << k << ") : " << Oefening3a(x) << endl 

k++; 


} 


Dit geeft als output: 


f(x) met x=10 A (-0): 0.158529015192103500000000000000 
f(x) met x=10 A (-l): 0.1665833531718476600000000000000 
f(x) met x=10 A (-2): 0.1666658333353170900000000000000 
f(x) met x=10 A (-3): 0.166666658333286700000000000000 
f(x) met x=10 A (-4): 0.1666666665852405100000000000000 
f(x) met x=10 A (-5): 0.1666666669955220400000000000000 
f(x) met x=10 A (-6): 0.1666666570693547400000000000000 
f(x) met x=10 A (-7): 0.1666639687323639700000000000000 
f(x) met x=10 A (-8): 0.1664054747909410900000000000000 



f(x) met x=10 A (-9): 0.2019483917365789700000000000000 
f(x) met x=10 A (-10) : 0.1415461897416452900000000000000 
f(x) met x=10 A (-ll): 0.1814452431681659900000000000000 
f(x) met x=10 A (-12): 0.1554369993887551900000000000000 
f(x) met x=10 A (-13): 0.1151016729396209700000000000000 
f(x) met x=10 A (-14): 0.1663955637223883700000000000000 
f(x) met x=10 A (-15): 0.1809989821226363200000000000000 


3b) 


De lim v ^ 0 ƒ (x) berekenen we als volgt. 


lim/(x) = lim 

x—>0 x—>0 


r x- sin x'' 


V 


) 


Hiervan de limiet nemen is echter niet mogelijk, we hebben een regel/eigenschap nodig om 
dit te verkrijgen. De regel hiervoor nodig is de regel van L’Hopital: 

r. 


lim| 

*->() 


r x- sinx A 


V 


) 


= lim| 

x->0 


(x - sin x) 

w 


We gaan dit dus afleiden. We verkrijgen: 
O-cosx^ 


lim| 

x->0 


V 


3x 


) 


We passen opnieuw de regel van L’Hopital toe: 


lim| 

x-»0 


(l - cos x) 


= lim[ 

x —>0 


^sin x ^ 


V 6x j 


We passen nog een laatste keer de regel van L’Hopital toe: 
(sin x) 


lim 

x->0 


(6x) 


liml 

x->0 


^cosx^ 


V o j 


Hiervan de limiet berekend is gelijk aan: 


lim/(x) = lim 

x—>0 x—>0 


^cosx^ 


V O J 


- = 0,16666666666666666666666666666667 
6 



3c) 


De functie is als volgt: 


r . . v - sin x 

ƒ (*) = - 3 - 

X 

Waarbij x heel erg klein wordt, namelijk machten van 10 k metk = 1,2, ... , 15. 

Om te laten zien hoe x verloopt: 

1 

0,1 

0,01 

0,001 

0,0001 

0,00001 

0,000001 

0,0000001 

0,00000001 

0,000000001 

0,0000000001 

0,00000000001 

0,000000000001 

0,0000000000001 

0,00000000000001 

0,000000000000001 

We dienen onze berekening zodanig te organiseren dat we rekening houden met een aantal 
basisregels. Zo zien we dat x 3 minder klein is dan sin x voor hele kleine waarden van x. Dit 
kunnen we gebruiken. Bovendien is getallen aftre kk en steeds zeer gevaarlijk, in het bijzonder 
indien het een relatief groot getal aftrekken is van een relatief klein getal. Ook een relatief 
groot getal delen door een relatief klein getal is niet goed en beter te vermijden. 

Bij het aftre kk en moeten we rekening houden dat het kleinste getal afgetrokken wordt van het 
grootste getal indien mogelijk. 

Het efficiëntste om het uit te rekenen is dus (in C++ code): 

doublé Oefening3 (doublé x) { 

return ((x + sin(-x))/pow(x, 3)); 

} 

Minder efficiënt zijn dus: 

doublé Oefening3b (doublé x) { 

return (1/pow(x,3))*(x-sin(x)); 

} 

doublé Oefening3b (doublé x) { 

return ((1/pow(x,2)) - (sin(x)/pow(x,3))); 

} 

doublé Oefening3b (doublé x) { 

return (1-(sin(x)/x))*(1/pow (x,2)) ; 

} 



doublé 0efening3b (doublé x) { 

return (1/pow(x,3))*(sin(-x)+x); 

} 

doublé Oefening3b (doublé x) { 

return ((sin(-x)/pow(x,3)) + (1/pow(x, 2))) ; 

} 

doublé Oefening3b (doublé x) { 

return ((sin(-x) + x)/pow(x, 3)); 

} 

doublé Oefening3b (doublé x) { 

return ((sqrt((1-cos(2x)) /2) + x)/pow(x, 3)); 

} 


Om tot de juiste oplossing te komen dienen we de sinus om te vormen naar de 
reeksontwikkeling van de sinus. Hierdoor krijgen we betere resultaten uiteindelijk en krijgen 
we minder fouten waardoor ons resultaat wél goed wordt. 

De reeksontwikkeling van de sinus is als volgt: 

oo 2 «+l 

sin* = V(-1)" —- voor alle x 

6? (2/i + l)! 


We dienen hiervoor in C++ ook nog de faculteit te programmeren, dit is eenvoudig te doen 
met recursie (maar kan ook indien nodig zonder recursie gemaakt worden en dan met een lus). 



Dit in C++ omgezet geeft: 

doublé faculteit (int i) { 
if (i <= 0) { 

return 1; 

} 

else 

{ 

return i * faculteit (i-1); 

} 

} 


doublé SinReeks (doublé x) { 
int i; 

doublé oplossing; 
doublé tussenoplossing; 
doublé teken; 

teken = 1; 

i = 1; 

oplossing = 0; 
while (i < 9) { 

tussenoplossing = pow(x, i); 

tussenoplossing = (tussenoplossing / faculteit(i)) 
if (teken == 2) { 

teken = 1; 

tussenoplossing = - tussenoplossing; 

} 

else 

{ 

teken = 2; 

} 

oplossing = oplossing + tussenoplossing; 
i = i + 2 ; 

} 

return oplossing; 

} 


doublé Oefening3beste (doublé x) { 

return ((SinReeks(-x) + x)/pow(x, 3)); 

} 


int main () { 


while (k <= 15) { 

x = (1 / pow(10, k)); 

cout << "f(x) met x=10 A (-" << k << ") : " << 

Oefening3beste(x) << endl; 
k++; 

} 

} 


En dit geeft de meest correcte oplossing. 



